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Abstract. We obtain exact results for the acceptance ratio and mean squared 
displacement in Monte Carlo simulations of the simple harmonic oscillator in D 
dimensions. When the trial displacement is made uniformly in the radius, we 
demonstrate that the results are independent of the dimensionality of the space. We 
also study the dynamics of the process via a spectral analysis and we obtain an accurate 
description for the relaxation time. 
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1. Introduction 

Since the original Metropolis algorithm appeared five decades ago, countless studies 
have employed the technique to evaluate the thermodynamic properties of model 
systems lMen and Tildesley, 19871 [Binder, 1997| [Frenkel and Smit, 20021 . The essence 



of the method is to generate a sequence of configurations that represent a given 
thermodynamic ensemble, often the canonical ensemble. Properties of interest are then 
obtained as averages over the configurations. At each step of the simulation a trial 
configuration is obtained from the current one by making a random displacement in 
the configuration space. This might correspond to, for example, displacing a randomly 
selected particle. The trial configuration is either accepted or rejected with a probability 
given by the appropriate Boltzmann factor for the ensemble. In case of rejection, the 
current configuration is retained for use in evaluating the properties of interest. 

Since many applications of MC are computationally intensive, a much addressed 
issue has been the optimization of the simulation with respect to one or more control 
parameters so that the configuration space is sampled in the most efficient way. 
Bouzida, Kumar and Swendsen (BKS) [Bouzida et al., 1992[ Swendsen, 2002 , as part 



of a program aimed at improving the efficiency of MC simulations of biomolecules, 
performed numerical studies of the simple harmonic oscillator (SHO) where the 
convergence of the simulation depends on the maximum displacement. There is no 
unique measure of efficiency, but two simple choices are the mean squared and mean 
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absolute displacements. As the maximum displacement tends to zero or infinity it is 
clear that the average value of both of these quantities tends to zero since, in the first 
case the particle does not move, while in the second all attempted moves are rejected. 
Thus both quantities have a maximum for some intermediate value of the maximum 
displacement. 

It is also useful to consider the dynamical process associated with an MC simulation, 
even though it does not correspond to the actual dynamics of the system. One can, 
for example, calculate various time correlation functions that can be used to develop 
alternative efficiency criteria |Kolafa, 1988HMountain and Thirumalai, 1994| . 

BKS |Bouzida et al., 1992| performed numerical studies of the SHO in one, two and 
three dimensions and examined the acceptance ratio Pace (i-e., the fraction of accepted 
trial configurations), mean squared, < (Ax)^) >, and mean absolute, < |Ax| >, 
displacements as a function of the maximum displacement, 6. They found that the 
acceptance ratio decreases approximately exponentially for small to intermediate values 
of S and then inversely for larger values. In one dimension they found that the maxima 
in < (Ax)^ > and < |Aa;| > occur at Pace = 0.42 and Pace = 0.56, respectively. In 
higher dimensions the results depend on how the jump is made. BKS considered two 
cases: in one the jumps are performed uniformly to any point in a spherical volume of 
radius S centered on the current position, a choice that favors larger radial displacements 
for dimensions D > 1. In a second method, the jumps were sampled uniformly in the 
radius (and randomly in the orientation) so that all radial displacements are equally 
probable. In the former case BKS observed that, for a given 6, Pace decreases as a 
function of D, while for uniform radius sampling the numerical results suggested that 
Pace is independent of D. 

In addition to these static properties, BKS also examined the correlation time, r, 
of the energy-energy correlation function. They observed a minimum correlation time 
for an acceptance ratio of approximately 50%. 

Here we present an analytical study of the SHO in arbitrary dimension. We obtain 
exact expressions for the acceptance ratio and the mean squared and mean absolute 
displacements as functions of the maximum displacement 6. We show that when the trial 
jump is selected uniformly in the radius, the results are independent of the dimension. 
We also present an analysis of the dynamics of the process. 



2. ONE DIMENSION 



We first investigate the case of a SHO in one dimension, whose potential energy is given 
by V{x) = kx'^/2 where x is the position and k is the stiffness constant. 

In a standard Metropolis Monte Carlo Simulation one makes a trial move with 
a uniform random displacement selected between —S and 6. The dynamical process 
generated by the successive trial moves of the Monte Carlo simulation can be written 
dP{x,t) 



dt 



1 r I r 

-— dhW{x ^ X + h)P{x,t) + — dhW{x + h^ x)P{x + h,t) {1] 
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^5 for a harmonic oscillator with 
uniform volmne sampling in D = 1,2, and 3 dimensions (full, dashed, and dotted 



Figure 1. Acceptance ratio Pace versus ^ — \/ 2 

IK 

lines respectively). For a uniform radius sampling, all curves coincide with the one- 
dimensional result (full line). 



where W{x —>■ x + h) denotes the transition rate from the state x to the state x + h 
and P{x, t) is the probabihty to find the oscillator at position x at time t. To ensure 
the convergence towards equilibrium, a sufficient condition is given by detailed balance, 
which is expressed as 



W{x x + h) _ Peq{x + h) 

W{x + h X)~ Peq{x) 



where 



cexp{—j3V{x)) 



-00 dxPgqi^X^ 



(2) 

(3) 
1. One 



and c = ^pk/2iT is a normalization constant ensuring that /! 
solution of Eq. ^ is the Metropolis rule, 

W{x x + h) = min(l, exp(-/3(\/(x + h) - V{x))). (4) 

Most of the properties of the ID SHO can be obtained analytically. For instance, 
the acceptance ratio, which is the number of accepted trials over the total number of 
trials can be expressed as 



PacciS) 



-00 2o 



rS 

/ dhW{x ^ X + h). 

J —S 



(5) 
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Figure 2. Mean squared displacement < (Ar)'^ > ^ versus the acceptance ratio Pace 
for uniform volume sampling in one, two and three dimensions (full curve, dashed, and 
dotted lines, respectively) . For a uniform radius sampling all curves coincide with the 
one-dimensional result (full line). 



In this equation exp{—i3V{x))dx is the probabihty that the oscillator is between x and 
X + dx, dh/26 is the probability of selecting a random displacement between h and 
h + dh and W{x ^ x + h) is the probability of accepting the trial displacement (given 
by Eq. (@))). Integration over the allowed values of x and h then gives the average 
acceptance probability. Since displacements to the left and right are symmetric, we 
need consider only one direction. For displacements to the right, Eq. (0), can be written 
as 



d6 



For X < -5/2, W{x 
One thus obtains 

di 



dxe~^^^''^W{x ^ x + 5). 
x + 5) = 1 and for x > -S/2, W{x - 



(6) 



x + S) 



-f3k{{x+S)^-x^)/2 



1 — erf 



(7) 



where erf{x) is the error function and ^ 



2 



Using the initial condition, i.e. 



-Pacc(O) = 1, the solution of the differential equation ((7j) is 



(8) 



The function is plotted in Fig. ^ 
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Figure 3. Mean absolute displacement < | Ar| > y ^ versus the acceptance ratio Pace 
for a uniform volume sampling in one, two and three dimensions (full curve, dashed, 
and dotted lines, respectively). For a uniform radius sampling all curves coincide with 
the one-dimensional result (full line). 



The mean squared displacement, < (Ax)^ >, and the mean absolute displacement, 
< I Ax I >, defined as 

dxexp{-f3V{x))— dhh^W{x ^x + h), (9) 

-co 2d J —5 

/+00 \ p5 

dxexp{-(3V{x))— dh\h\W{x ^ X + h). (10) 
-oo 2o J -5 

can be quite simply obtained from the generating function 

Z{X) = c dxexp{-f3V{x))— dhexp{-X\h\)W{x ^ x + h), (11) 



26 J-5 

by the derivatives with respect to A and set A = 0, i.e., < |Ax| >= —{dZ{\)/d\)x=o, 
< {Axf >= {d'^Z{X)/dX^)x=o (and, of course. Pace = Z{X = 0)). 

By multiplying both sides of Eq. (|TT| by ^ and next differentiating with respect to 
^, one obtains 

««^(«))-exp(-,/lA«Ul-e./m) (12) 



V V /^^ A V2 



which, after defining A = y^A, leads to 



Z(A,0 = i 



1 - exp(-AO (l - erf ( ^ j j - exp(P) Lrf ( f + A J - er/(A) 



(13) 
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By taking the derivatives of the above formula with respect to A and evaluating 
the resulting expressions at A = 0, it is easy to show that the mean squared and mean 
absolute displacements are given by 




< lAxI > 




e 



erf (I) 





exp 



(14) 



(15) 



The maximum in the mean squared and mean absolute displacements occur for 
^ = 2.61648 and ^ = 1.76332, which corresponds to acceptance ratio values of 
Pace = 0.41767 and Pace = 0.558239 in agreement with the numerical results of 
BKS IBouzida et al., 19921 : see Figs. El and El 



3. D dimensions 



We show here that the acceptance ratio, mean squared displacement and other quantities 
of interest can be obtained exactly in any dimension. Note that the term "volume" 
should be interpreted as the hypervolume in D dimensions, e.g., area in 2D and volume 
in 3D. For simplicity, we derive exact expressions in odd dimensions, but similar results 
can be obtained in even dimensions. 



3.1. Uniform volume sampling 

In D dimensions, the acceptance ratio is expressed as 

PaceiS) = f d^rexp(-/5A;rV2) / rf^hmin(l, exp(-/5A;(|r + hp - r2)/2)) (16) 

O^Vd J J\h\<5 

where Vd = 7r^/^/r(D/2 + 1) is volume of the sphere of unit radius in D dimensions 
and 



CD 



-I /,3/.\^/2 



DVd drr^-^e-^'^'/^] =(gj . (17) 



In odd dimensions, the derivative of Eq. (fTHIl with respect to 5 can be written 
explicitly by using generalized spherical coordinates 

*^!^ = ££y<i«rexp(-/JfcrV2) 

s"-^ j da min(l,exp(— /5fc(r5cos(/)i + 5^/2))), (18) 

where dQ = {Ufji{sH<Pj))^'^~^d(f)j)d(j)D-i such that J dQ = DVd- The first D - 2 
variables (pj are integrated from to vr, whereas (pD-i is integrated from to 2tt. If 
we denote u = cos(0i), perform integration over (p2 ■ ■ ■ (pD-i, and introduce the variable 
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V = r\J^ Eq. (fTH|) can be rewritten as 

r(f)/_^lC/M(l-M2)(^5-3)/2 

X v^'-^dve-''' J ^du{l-u^)'^''-''^/^mm{l,exp{-{2^vu + e))) (19) 

where min(l, exp(— (2^wm + ^^))) = exp{{—{2^vu + ^^)) for v < ^/2 with — 1 < m < 1 
and for v > ^/2 with -^/{2v) <u < 1 and min(l, exp(-(2^WM + ^^))) = 1 for r > ^/2 
with -1 < M < -^/{2v). Using that /\ c/m(1 - m2)(^-3)/2 ^ £iV)^^ ^q. ((121) then 
becomes 



dvv'^-'e-"" J ^ du{l - m2)(«-3)/2 exp(-(e2 + 2^vu)) 



r+oo 

+ / v'^-'dve-''^ 

Ji/2 



■ 1 



+ / du{l-u'f'^ex^{{-{e + '^ivu))) 

J-(/i2v) 

After some calculation (see Appendix A) one obtains. 



(20) 



which gives, for instance, in three dimensions 

PUO -I- erf (f ) + -A_ (l - (l 4) exp (4)) (22, 

Figure H shows the acceptance ratio Pace versus 6 in one, two and three dimensions. 
A similar calculation for generating function Zd{X,^) leads to 

'J^^Mm.ne-^.-^^[,-.rf(f)). (23) 

and to 

ZdCKO = D{-0'~^£^^Zn=,{\0. (24) 

where the expression for Z£)=i(A,^) is given in Eq. (fT^ . Since Pacc(0 = Z£,{X = 0,^), 
it follows from Eq. that the acceptance ratio in D dimensions is equal, up a 
factor D(-0^~^yf , to the mean squared displacement < (Ax)^ > in one dimension: 
compare Eq. ^ to Eq. (fi3|l . 

Although straightforward, the algebra rapidly becomes tedious, and we only 
illustrate the results by giving the expression of the mean squared displacement 
< (Ar)^ > in three dimensions 

f <(Ar)^> =^^z.^3(A,OIa.o 
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The mean squared and mean absolute displacement are plotted versus the acceptance 
ratio Pace in one, two and three dimensions in Figs. |21and El Note that the maximum 
is shifted to the left, i.e., to the smallest values of the acceptance ratio, when the space 
dimension increases. 



3.2. Uniform radius sampling 

The acceptance ratio Pacc,w in D dimensions can be expressed as 
Pacc,M=CDf c/^rexp(-/5A;rV2) 

JD 

I d^hP^(/i)min(l,exp(-/3A;(|r + hP-r2)/2)) (26) 

J\h\<5 

where Pw{h) is the weighted probability. For a uniform distribution in radius, 
h^~^Pu,{h) = {DVdS)~^. Using the method developed in the above section, it is 
straightforward to obtain that 

^^M^^(l-.r/(|)) (27) 

which shows that the acceptance ratio is the same whatever the dimension and explains 
the data collapse observed in [Bouzida et al., 1992| . Similarly, the generating function 
Z(A,^) can be shown to obey to the differential equation 

independently of the dimension D. This proves that Z{X,C,) and all moments such 
as < |Ar| > and < (Ar)^ > are independent of dimension, as numerically found by 
BKS IBouzida et al., 19921 



4. Dynamic behavior 

In addition to the exact results for the static properties presented above, we have 
investigated the dynamic behavior of the SHO by using numerical and analytical 
approaches. For simplicity, we discuss only the unidimensional case, but the approach 
can be generalized to D dimensions. 

The master equation describing the dynamical evolution of the system during the 
Monte Carlo simulation, Eq. (^, can be formally written 

where L is a linear operator acting on P and the Metropolis rule, Eq (ji}, is used for the 
transition rate. We consider the spectrum of eigenvalues A of L. Denoting P\{x) the 



Optimum Monte Carlo Simulations: Some Exact Results 



9 



eigenfunction associated with A and introducing fx{x) via Px{x) = Peq{.x)f\{x), where 
Peq{x) is given in Eq. Q, one can express the eigenvalue equation 

XPeq{x)fx{x) = L{P,,{x)fx{x)) (30) 

as 

A/a(x) = ^ I' rf/iMzn(l,e-'3(^(^+'^)-^(^)))(/A(a;) - + h)). (31) 

Multiplying both sides by exp{—pV{x))f^{x), where the star denotes a complex 
conjugate, and integrating over x then gives 

^ 1 /+^rfx/^,rf/iMm(e-/^^(-),e-/^^(-+^))|A(x + /i) - Mx)\' 

45 i:!:^dxe-pn^)\MxW ' ^ ' 



As anticipated for a Markov process satisfying detailed balance, one deduces from 
the above formula that all eigenvalues are real and positive; the smallest eigenvalue is 
Ao = and it is associated with fo{x) = constant ^ 0. One need consider only real 
eigenfunctions. Moreover, the eigenvalues can be sorted according to the symmetry of 
the associated eigenfunctions: it is easy to check that the eigenfunctions are either even 
or odd functions of x, due to the fact that the potential V{x) is an even function of x. 

Any solution of the master equation can be expanded as 

P{x,t) = P,,{x){l + J2cxfx{x)e-^'), (33) 

A>0 

and a similar expansion applies to the conditional probability P{x,t\xo,0) from which 
one can compute any time-dependent correlation function. The long-time kinetics 
governing the approach to equilibrium in P{x, t) and in any correlation function is 
characterized by the smallest non-zero eigenvalue for which the amplitude, i.e., the 
projection of P(x, 0), or of the dynamic observable, onto the relevant eigenfunction, 
does not vanish. 

Since, according to Eq. the eigenvalues are expressed as the ratio of 

two positive quadratic functionals (that in the denominator being also definite), 
one can use the Rayleigh-Ritz procedure to find a variational upper-bound for 
the eigenvalues [Dettman, 19621 [Arfken, 1985| . Consider first the smallest non-zero 
eigenvalue Ai. For any real function (j){x) which is both normalized and orthogonal 
to fo{x), i.e., satisfies 



oo 

2, 



Peq{x)(j){xydx = 1 (34) 

Pegix)(p{x)dx =0, (35) 
one has the inequality 



oo 
oo 



oo 



1 



hoo 



<5 



Ai < Ai[</)] = ^ / dx c//iMm(e-'^^("),e-^^("+'^))|</)(x + /i) -</.(a;)|l (36) 

4o J-oo J-S 

A convenient choice of trial functions is provided by linear combinations of Hermite 
polynomials H^i^) (where, as in the previous sections ^ = J ^6), with n > 0, since 
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these form, up to a trivial multiplicative factors, an orthonormal basis with respect to 
the weight function exp{—j3V{x)) with V{x) = {l/2)kx^. For Ai, which is associated 
with an odd eigenfunction, one need consider only the odd polynomials if2n+i(0? n >0. 
The simplest estimate of Ai is provided by taking 

m = 4^ = (37) 



12^ 

which gives 




With this choice of 0(0) ^il^P] simply reduces to the mean squared displacement 
< (Ax)^ > multiplied by (^^j (see Eq. (fT^V One then derives from section 2 that 
Xi[4>] versus ^ passes though a maximum for ^ ~ 2.611648, which corresponds to an 
acceptance ratio of Pace = 0.41767. 

A better estimate of Ai can be obtained by using a linear combination of Hi{^) and 

^3(0: 



(39) 



'2^ ^48^ 

where only one independent parameter 6 appears due to the normalization condition. 
The best bound is determined by minimizing the expression Ai[0] with respect to 6: 

= 0, (40) 

The result is a lengthly algebraic formula that is plotted in Figs. |3] and together 
with the expression in Eq. ()38|) . 

An improved estimate of Ai can be derived by noting that at large ^, Ai is inversely 
proportional to ^. Actually, one can show that this is true for all eigenvalues except 
Ao = 0. By considering Eq. (p?T| in the limit where x goes to zero, one arrives at the 
result (see Appendix B) 



A(0~^ + O(e-^^ (41) 

valid for large ^. Note that the correction terms are very small as soon as ^ > 3. One 
can build an estimate of Ai(0 by using the piecewise function that is equal to ^ for 
^ > ^* and is equal to Ai [0] obtained for a linear combination of }i\ and (see above) 
for < where ^* is the value at which Ai[0] = A(0 is then maximum for 
^ = ^* ~ 2.35; the corresponding value of the acceptance ratio is Pace — 0.56. The 
estimate is shown in Figs. |3] and 

In order to compare our results to the BKS paper [Bouzida et al., 1992| , it is 
necessary to calculate the second eigenvalue A2. The trial function is then chosen 
in a subspace orthogonal not only to /o(0 = constant but also to the eigenfunction 
associated with Ai. A convenient choice is provided by (normalized) linear combinations 
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Figure 4. Ai versus ^ . The full curve was obtained by numerieal diagonalization of the 
master equation. The dashed curve corresponds to the zeroth-order estimate, Eq. I|38(l . 
the dotted curve corresponds to the solution of the first-order trial function, Eq. l|4Q|l l. 
and the dash-dot curve corresponds to the exact asymptotic behavior, Eq. H41|l . 



of the even Hermite poljTiomials H2n{C) with n>l. An estimate of A2 can be obtained 
by using a hnear combination of -^2(0 ^"^^ Hi{^): 



sm{9) 



(42) 



_2^2v^ 8^60F 
and the resuh is shown in Figs. IHl and [71 together with the zeroth-order approximation 
obtained with only -^2(0 ^^e improved estimate taking into account the large-^ 
behavior. 

In addition to the above analytical estimates, we have also performed a numerical 
study of the spectrum of eigenvalues of the master equation, Eq. ((H). The latter has 
been discretized in x-space by taking a constant step size A, which leads to a matrix 
form, 



dP{t) 
dt 



-WP(t) 



(43) 



where P is a vector with components Pi{t) = P{xi, t) and the elements Wij of the matrix 
W are such that 

i+Nh 
j=i~Nh,j^i 
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2 



1.5 



0.5 



"0 1 2 3 4 5 

Figure 5. Same as Figure ^ except versus ^. 

W,, = -AW{j ^ z) (44) 

and W{i j) = Mm(l, e"^*-^*-"^^-'"^*-^'-*''). The eigenvalues A can then be obtained 
via an exact numerical diagonalization of the matrix W. In practice, convergence is 
obtained for a unidimensional lattice of 400 sites, where the x-range is [—10, 10], and ^ 
goes from to 5. One checks that the lowest eigenvalue is equal to zero and corresponds 
to the equilibrium state, and that all other eigenvalues are real and strictly positive and 
behave as \/tt/{2^) for large enough ^. The results for the first non-zero eigenvalue Ai 
and A2 are displayed in Figs. Eland IHl respectively. One can see that the best analytical 
estimate described above (linear combination of two Hermite polynomials plus exact 
asymptotic behavior at large ^), is in excellent agreement with the numerical value in 
both cases. 

The above analysis allows us to derive by analytical means the simulation result 
obtained by BKS |Bouzida et al., 1992| for the dependence of the characteristic time r of 
the energy-energy correlation function on the acceptance ratio: approximating r by I/A2 
(since the energy V{x) is an even function of x, its projection on the first eigenfunction 
associated with Ai vanishes), using for A2 our best analytical estimate, and combining 
this with the exact result for the acceptance ratio Pace in section 2 lead to the full curve 
plotted in Fig. |H1 the time r is minimum for the acceptance ratio close to 0.47, as found 
by BKS IBouzida et al., 1992| (~ 0.50). 

In Fig. ISl we have also plotted 1/Ai versus the acceptance ratio: it is minimum 
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Figure 6. A2 versus ^ . The full curve was obtained by numerical diagonalization 
of the master equation. The dashed curve corresponds to the zeroth order estimate, 
{6 = in Eq. (|42(l ). the dotted curve correspond to the solution of the first order 
trial function, Eq. 142|l and the dash-dot curve corresponds to the exact asymptotic 
behavior, Eq. iglj. 



for Pace = 0.45. Note that using the results of this and the preceding sections, one 
can rigorously show that the correlation time r, no matter how it is precisely defined, 
diverges as 1/Pacc when Pace and as 1/ (1 - Pace) when Pace 1- 

5. Conclusion 



We have obtained exact results concerning Metropolis Algorithms for the displacement 
of a particle in the simple harmonic potential. Our analysis provides a theoretical 
explanation of the numerical results obtained by BKS |Bouzida et al., 1992| . In 
particular, we show that the results become independent of the space dimension when 
the successive trial moves are sampled according to a Metropolis algorithm with a 
uniform distribution in radius (instead of volume). This rationalizes the search for 
efficient Monte Carlo methods for the simulation of systems with intrinsic inhomogeneity 
and anisotropy such as biological molecules |Bouzida et al., 1992 Swendsen, 2002| 
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Figure 8. Inverse of the first nonzero eigenvalues, f/Ai (dotted curve) and I/X2 (full 
curve), versus the acceptance probability Pace 
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Appendix A. Acceptance probability 

The derivation of Eq. (pT|l can be done from Eq. piU) after some manipulations whose 
details are given here. Let us denote Id and Jo as 

Id = dw^'-h-'"' J ^ du{l - m2)(^-3)/2 exp(-(e' + 2^vu)) (A.l) 



+00 



-?/(2f) 

1 



Jd = / dvv^-^e-'" 

du{l - exp((-(e' + 2ivu))) 



(A.2) 



Changing the variable u to y = vu and v to t = v"^ — y"^ leads to the following relations 
Jo 2 Jo 



J5/2 2 Jo 

Jo 2 Jei^-y'^ 

Jo 2 Jf 



-00 ^^rv^/ /n, I C\2'l 



+ / ■ rf/"P(-(^ + 0^) dt t(--)/^e- (A.4) 



/C/2 2 Jo 

Using that J+°^ dtt^^-^^/^e'^ = T{{D - l)/2), one obtains that 



2 \ 2 J \ V2. . 

Inserting Eqs. (fX;5|l in Eq. ^ leads to Eq. (j2T]l . 

Appendix B. Asymptotic behavior of the eigenvalues for large 6 

Consider Eq. (jHlj) in the limit x ^ ; one obtains after rearranging the various terms: 



1 ...1 , .^2 2^ 



+ - / dh{l - exp(^((a; + /i)^ - x'mU^) - fxix + h)). (B.l) 
2d J-2x 2 

The second term of the right-hand side of Eq. ()B.1|) is at most of order x^\fx{x)\ and is 
always negligible so that one can rewrite Eq. ()B.1|) as 



X6 - dhe^''" + Oix")] fx{x) ~ dhexpi^ix + hf)fx{x + h).{B.2) 
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Shifting the variable from h to x + h in the integral of the r.h.s of Eq. ()B.2|) and 
using the orthogonality of f\{x) to /o(a;) = constant for all non-zero eigenvalues A, i.e., 
dx exp{^^x'^) f\{x) = 0, leads to the following expression, 



\6 



+ 0(x2) /a(x) ~ - 



dhe-^h{h)+ / dhe-^h{h) 

x+5 J ~oo 



for any non-zero A. The r.h.s. of Eq. ()B.3|) can be Taylor expanded, which gives 



(B.3) 



\6 



TC 



1 



+00 



-I3kh 



dhe^{h{h) + M-h)) 



-I3kh 



xe^iMh)-M-h)) + 0{x') 



(B.4) 



If the eigenfunctions fx is an even function of x, one then derives that f\{x 
/a(0) + O(x2) with /a(0) ^ and 

-Pkh\fx{h) 



X6 



+00 

which after introducing ^ = \J^5 can be rewritten as 



(B.5) 



A 



(B.6) 



If the eigenfunction is an odd function of x, one has that fx{x) = /a(0)x(1 + 0{x^) with 
MQ) ^ and 

(B.7) 



A = — — exp(-^ J 



2^ /a(0) 

From Eqs. ()B.6|) and ()B.7|) . one immediately obtains that all non-zero eigenvalues behave 



as 



A~^ + 0(exp(-a) 
when ^ +00, since f\{x) diverges more slowly than e^^ when x — > -|-oo. 



(B.8) 
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